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In  studies  of  radiation  from  partially  coherent  sources  the  so-called  quasi-homogeneous  (QH)  model  sources 
have  been  very  useful,  for  instance  in  elucidating  the  behavior  of  fields  produced  by  thermal  sources.  The 
analysis  of  the  fields  generated  by  such  sources  has,  however,  been  largely  carried  out  in  the  framework  of 
scalar  wave  theory.  In  this  paper  we  generalize  the  concept  of  the  QH  source  to  the  domain  of  the  electromag¬ 
netic  theory,  and  we  derive  expressions  for  the  elements  of  the  cross-spectral  density  matrix,  for  the  spectral 
density,  the  spectral  degree  of  coherence,  the  degree  of  polarization,  and  the  Stokes  parameters  of  the  far  field 
generated  by  planar  QH  sources  of  uniform  states  of  polarization.  We  then  derive  reciprocity  relations  analo¬ 
gous  to  those  familiar  in  connection  with  the  QH  scalar  sources.  We  illustrate  the  results  by  determining  the 
properties  of  the  far  field  produced  by  transmission  of  an  electromagnetic  beam  through  a  system  of  spatial 


light  modulators.  ©  2005  Optical  Society  of  America 
OCIS  codes:  030.1640,  260.2110,  260.5430. 

1.  INTRODUCTION 

The  so-called  quasi-homogeneous  (QH)  sources  (see  Ref. 

1,  Sec.  5.3.2)  are  important  models  for  many  partially  co¬ 
herent  sources  found  in  nature  or  developed  in  laborato¬ 
ries.  Use  of  such  models  has  clarified  the  behavior  of  Lam¬ 
bertian  sources2  as  well  as  of  secondary  sources  that  are 
encountered  in  scattering  theory. 3-7 

Up  to  now  almost  all  investigations  concerning  QH 
sources  and  the  fields  that  they  generate  have  been  based 
on  scalar  wave  theory.  Perhaps  the  most  important  and 
best-known  result  derived  from  the  scalar  QH  model  is 
the  van  Cittert-Zernike  theorem  and  its  generalization 
(see  Refs.  1  and  8),  which  give  a  simple  expression  for  the 
degree  of  coherence  at  a  pair  of  points  in  the  field  propa¬ 
gated  from  such  a  source.  When  applied  to  points  in  the 
far  zone  the  generalized  van  Cittert-Zernike  theorem  ex¬ 
presses  a  reciprocity  between  the  intensity  distribution 
across  the  source  and  the  degree  of  coherence  of  the  far 
field.  There  is  a  second  less  well-known  reciprocity  rela¬ 
tion  for  fields  generated  by  such  sources  that  expresses 
the  spectral  density  distribution  in  the  far  zone  in  terms 

1084-7529/05/112547-10/$15.00 


of  the  spectral  degree  of  coherence  of  the  source  (see  Ref. 

1,  Sec.  5.3.2). 

In  this  paper  we  generalize  the  concept  of  the  quasi- 
homogeneous  source  to  the  domain  of  the  electromagnetic 
theory  with  the  help  of  the  recently  introduced  2X2 
cross-spectral  density  matrix  of  the  electric  field. 9-11  We 
then  show  that  a  beam  generated  by  an  electromagnetic 
QH  source  that  is  uniformly  polarized,  i.e.,  that  has  the 
same  state  of  polarization  at  each  point,  obeys  two  reci¬ 
procity  relations  analogous  to  those  that  are  familiar  for 
scalar  QH  sources.  We  also  derive  expressions  for  the 
spectral  degree  of  polarization  and  the  spectral  Stokes  pa¬ 
rameters  of  the  far  field  that  are  generated  by  such 
sources.  We  illustrate  the  results  by  plots  that  show  the 
behavior  of  the  far  field  produced  by  transmission  of  such 
a  beam  through  a  system  involving  spatial  light  modula¬ 
tors. 

2.  QUASI-HOMOGENEOUS 
ELECTROMAGNETIC  SOURCE 

Consider  a  fluctuating,  planar,  secondary  electromagnetic 
source,  located  in  the  plane  z  =  0  and  radiating  into  the 

©  2005  Optical  Society  of  America 
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Fig.  2.  Illustration  of  the  concept  of  a  QH  source. 

half-space  z>  0.  We  assume  that  the  radiated  field  is 
beamlike,  propagating  close  to  the  z  axis,  and  that  the 
source  fluctuations  are  represented  by  a  statistical  en¬ 
semble  that  is  stationary,  at  least  in  the  wide  sense.  The 
second-order  correlation  properties  of  the  source  may  be 
characterized  by  the  2X2  electric  cross-spectral  density 
matrix:9 

W<0)(Pl  ,P2,®)  = 


^(puPiM  W<>1,P2,0>) 

W(y°J(Pl,p2,co)  WfJ(Pl,p2,co) 


(2.1) 


In  this  formula  Pi  and  p2  are  the  two-dimensional  posi¬ 
tion  vectors  of  two  points  Qi  and  Q2  in  the  source  plane 
(see  Fig.  1),  co  denotes  the  frequency,  and  the  elements  of 
the  W-matrix  are 

y^\pi,p2,M)  =  {E*{Pi,(ti)Ej{p2,(o))  (i i=x,y ;  j  =  x,y). 

(2.2) 

Here  E,  and  Ej  denote  the  Cartesian  components  of  a  typi¬ 
cal  member  of  the  statistical  ensemble  of  the  electric  field 


in  two  mutually  orthogonal  x  and  y  directions  perpendicu¬ 
lar  to  the  z-axis  (the  axis  of  the  beam),  and  the  angle 
brackets  denote  the  average,  taken  over  an  ensemble  of 
realizations  of  the  electric  field,  in  the  sense  of  coherence 
theory  in  the  space-frequency  domain  (Ref.  1,  Sec.  4.7.1). 

The  electric  cross-spectral  density  matrix  representing 
an  electromagnetic  QH  source  may  be  obtained  as  a 
straightforward  generalization  of  the  cross-spectral  den¬ 
sity  function  of  a  scalar  QH  source:1,8 


Wf\Pl,p2M  =  M0)((Pi  +  p2)/2,«)Vsj0)((p1  +  P2)/ 2,0.) 

X/4°’(P2-Pl>"), 

( i=x,y ;  j  =  x,y),  (2.3) 

where  Sf\P,co)=W(°){P,P,  co)  and  S^)(P,io)=W^\P,P,co) 

are  the  spectral  densities  of  the  components  E,  and  E .  and 
fo)  J 

p  (p2-pi,co)  represent  the  correlations  between  these 


components.12  It  is  assumed  that  S{0\p,co )  ( i=x,y )  varies 
much  more  slowly  with  p  than  p,)  .  (p' ,  co)  varies  with  p' 
(see  Fig.  2). 

The  spectral  density  »S(0)(p,  <o),  the  spectral  degree  of  co¬ 
herence  ?/()l(pi  ,p2,(u),  and  the  spectral  degree  of  polariza¬ 
tion  Pi0\p  ,  co)  of  an  electromagnetic  beam  are  given  by  the 
formulas9 

S(0\P,w)  =  Tr  W<0*(p,p,  co), 


V(0)(pl,P2,to) 


Tr  #<%>!, pa, «) 


Vs( 0)(p1,o«)Vs(0)(p2,o> ) 


P{0\p,co)  = 


4  Det  W{0\p,p,co) 
[Tr  W(0)(p,p,  a>)]2  ' 


(2.4) 


(2.5) 


(2.6) 


where  Tr  denotes  the  trace  and  Det  the  determinant.  For 
a  quasi-homogeneous  electromagnetic  source  whose  cross- 
spectral  density  matrix  W(0)(p1,p2,o>)  is  given  by  Eq.  (2.3) 
expressions  (2.4)— (2.6)  become 


S^(P,co)  =  S{°\P,co)+Siy°\P,co), 


(2.7) 


ym(pi,p2,M)  ■ 


Si0)((P!  +  p2)/2,oj)pS)(p2  -  Pi,  to)  +  S(°\(Pl  +  p2)/2 ,co)pfJ(p2  -  p1;  to) 
S(0)((Pi  +  P2)/2,<o)  ’ 


(2.8) 


V[Sf(p,to)  -  S‘0,(p,a>)]2  +  4Sf  (p,to)Sf  (p,co)|p^(0,t,)|2 
P{0\p,co)  = - — 


S(0)(p,to) 


(2.9) 


Formula  (2.7)  holds  generally;  for  a  QH  source, 
S^°Kp,co)  is  a  slow  function  of  p.  The  Stokes  parameters, 
which  are  useful  for  characterizing  the  state  of  polariza¬ 


tion  of  an  electromagnetic  field,  may  be  expressed  in 
terms  of  the  elements  of  the  cross-spectral  density  matrix 
by  the  formulas  (cf.  Ref.  13,  Sec.  10.9.3) 
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So0)(p,") 

=  W<i)(p,p,<U)  +  W^>(p,p,a»), 

(2.10a) 

Si °\p,co) 

=  WS)(p,P,®)-W^)(p,p,a.), 

(2.10b) 

s  $\p,co) 

=  2Re[W^(p,p,a,)J, 

(2.10c) 

S^0)(p,w) 

=  21m  [W?J(p,p,co)l 

(2.10d) 

Re  and  Im  denoting  the  real  and  the  imaginary  parts,  re¬ 
spectively.  It  will  be  useful  to  normalize  the  Stokes  pa¬ 
rameters  by  setting14 

s‘°W)  =  Sf  (p,<u)/S^0)(p,oj)  (j  =  1,2,3).  (2.11) 


The  normalized  Stokes  parameters  s*0>(p,oj)  (e  =  1,  2,  3)  as¬ 
sociated  with  the  QH  source  characterized  by  Eqs.  (2.3) 
are  readily  shown  to  be  given  by  the  formulas 


Sf(p,o))-Sf(p,o>) 

"Sf(p,a.)+S<0)(p,w)' 


(2.12a) 


s(20>(p,m)  = 


2  V'Sftp,  a,)  yj&Vjp,  o>)Re[pfflp,p,  o>)] 


(2.12b) 


S  30)(P,"): 


2  V  sf  (p,  o,)  ys<Q)(p,  o,)im[pg)(p,p,  q,)] 
Sf(p,6»)  +  Sf(p,a.) 


(2.12c) 


“y(")  =  < 


1  +  a(a i) 

\ja(w) 

1  +  a'(oj) 
<?(&)) 

1  +  a(a)) 


when  i=j  =  x 
when  i  +  j 
when  i=j  =y 


(2.16) 


The  spectral  degree  of  coherence  of  the  source  is  then 
given  by  the  expression 

v{0\pi  ,p2,(o)  =  axx(o)  )/“«(p2  -pi,to)  +  a  Ja>  )m1°v)(P2  -Pl,“)- 


(2.17) 


Expressions  for  the  normalized  Stokes  parameters  of  the 
source  are  derived  in  Appendix  A. 


3.  THE  SPECTRAL  DENSITY  AND  THE 
SPECTRAL  DEGREE  OF  COHERENCE  OF 
THE  FAR  FIELD.  TWO  RECIPROCITY 
RELATIONS 

It  can  be  shown  by  a  straightforward  analogy  of  the  argu¬ 
ment  used  in  deriving  Eqs.  (5.3-4)  and  (5.3-5)  of  Ref.  1 
that  the  elements  of  the  cross-spectral  density  matrix  of 
the  electric  field  in  the  far-zone  are  given  by  the  expres¬ 
sions  (see  Fig.  1) 


We  will  consider  only  sources  for  which  the  normalized 
spectral  Stokes  parameters  are  independent  of  p.  The 
spectral  degree  of  polarization  and  the  spectral  polariza¬ 
tion  ellipse  associated  with  the  polarized  portion  of  the 
beam  will  then  be  the  same  at  every  point  of  the  source. 
We  will  call  such  sources  uniformly  polarized.  In  Appen¬ 
dix  A  we  show  that  the  source  is  uniformly  polarized  at 
frequency  co  if  and  only  if  the  spectral  densities  S{°\p,(o ) 

/A}  X 

and  Sy  (p,io)  are  proportional  to  each  other,  i.e.,  if  and 
only  if15 

S™(p,o»)  =  a(o»)Sf(p,w)  (2.13a) 

and  if,  in  addition,  the  correlation  coefficient  p®y(p,p,oo) 
is  independent  of  position,  i.e., 

Px y(p,P,^)  =  Pxyiv).  (2.13b) 

It  follows  from  Eqs.  (2.4),  (2.13a),  and  (2.13b)  that 

Sf(p,m)  =  - - —  S(0W), 

1  +  a{u>) 


w<”VlSl  ,r2s2,&» )  =  (2tt&)2cos  Oi  cos  02 

xWi-k  S1L,kS2±,u) 


exp[ik(r2  -  rj)] 


r\r2 


where 


(3.1) 


ttf(fl’f2’“)  =  f2^  j  j  W»>(pi,P2,£ 


,<u)exp[-  i{f1  •  pj 
+  f2-p2)]d2pid2p2  (3.2) 


is  the  four-dimensional  spatial  Fourier  transform  of 
W^lpi  ,p2,oj).  In  Eq.  (3.1),  s,  and  s2  are  unit  vectors 
along  iq  and  r2  and  s1±  and  s2  are  the  projections,  con¬ 
sidered  as  four-dimensional  vectors,  of  the  unit  vectors  S! 
and  s2  onto  the  source  plane.  Each  integration  in  Eq.  (3.2) 
extends  over  the  source. 

On  substituting  from  Eq.  (2.15)  into  Eq.  (3.2)  we  obtain 
for  the  Fourier  transforms  of  the  elements  of  the  cross- 
spectral  density  matrix  of  a  uniformly  polarized  QH 
source  the  expressions 


a(co) 

S(y\p,w)  =  - - —Si0)(p,co).  (2.14) 

1  +  a((o ) 

On  substituting  from  Eqs.  (2.14)  into  Eq.  (2.3)  we  find 
that  the  elements  of  the  cross-spectral  density  matrix  of 
the  source  that  satisfy  conditions  (2.13a)  and  (2.13b)  are 
given  by 

W|f(Pl,P2,")  =  alJ(a))S{0}((p1  +  p2)/2,(o)p-j](p2  -  pi,w), 

(2.15) 


a.;(co)  f  f 

xPij  \p2  ~  Pi,  w)exp[-  i(f1  •  px 
+  f2  •  p2)]d2pid2p2 

=  +  f2,  o,)40)((f2  -  fj)/2, tu) 


where  the  coefficients  atJ( oj)  are 


( i=x,y ;  j  =  x,y),  (3.3) 
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where  a^ioj)  are  the  coefficients  defined  in  Eq.  (2.16)  and 


S(0)(f,  (o)  = 

40)(f>")  = 


— — ^  J  S(0)(p,<»)exp[-  ff-p]d2p, 

(27 r)2  J 

7Z  72  I  My0>(P,")exp(-if-p)d2p- 

(2ttJ  J 


(3.4) 


(3.5) 


Formula  (3.3)  is  derived  in  Appendix  B. 

On  substituting  from  Eq.  (3.3)  into  Eq.  (3.1)  we  obtain 
for  the  elements  of  the  cross-spectral  density  matrix  in 
the  far  zone  the  expressions 

Wl“'ViSi,r2S2,«) 

=  (2tt&)2COS  COS  #2 «y(a>) 

xS|0|Ws21-Sli),ffl) 


xAtij>)(^(s2x  +  Si±y2,&)) 


exp[l£(r2  -  rx)] 


r\r2 


(3.6) 


The  spectral  density  S^drs,  &j)  of  the  far  field  is  given 
by  the  expression 

S^\rs,(o)  ~  Tr  W^VrSjrSjiu)  =  (2ir£)2cos2  fiS(0)( 0,co) 


X\axx{io)fj,{x°J(ks±,(o)  +  ayy(<o)/l(°J(ks±,(o)]—, 

(3.7) 


positive  2  axis.  We  also  have,  on  taking  the  Fourier  trans¬ 
form  of  Eq.  (2.17),  that 

?7(0)(£sx,&>)  =  axx(a)flxx(ks  !  ,w)  +  ayy{u>)jfyy(k S±,a>). 

(3.8) 

On  substituting  from  Eq.  (3.8)  into  Eq.  (3.7)  we  obtain  for 
the  spectral  density  of  the  electric  field  in  the  far  zone  the 
expression 


S(co)(rs,o)  =  (2TTk)2cos28S{0)(0,co)yi0)(ks±,co)^:.  (3.9) 

rz 

Equation  (3.9)  expresses  a  reciprocity  relation  that  holds 
for  beams  generated  by  a  planar,  secondary,  uniformly  po¬ 
larized  electromagnetic  QH  source,  namely,  the  spectral 
density  of  the  far  field  is  proportional  to  the  product  of  the 
Fourier  transform  of  the  spectral  degree  of  coherence  of  the 
field  across  the  source  and  of  cos2  6/r2. 

Next  we  consider  the  spectral  degree  of  coherence  of 
the  far  field  at  points  r1=r1s1  and  r2=r2s2,  viz., 

77W(os1,r2s1,£o) 

Tr  #(”)(r1s1,r2s2,w) 

\  Tr  W^ViSj.rjSi,^)  VTr  WM(r2s2,r2s2,a>) 

(3.10) 


where  6  is  the  angle  that  the  unit  vector  s  makes  with  the 


With  the  help  of  Eq.  (3.6)  one  finds  that 


77W(riSi,r2s2,(o): 


axx(a})jl^x)(k(s.2j  +  s1±)/2,w)  +  ayy(co)p(yy)(k(s2±  +  s1±)/2,ai) 
'J‘xXx(w)gxx(ksl±’c°)  +  “xyMaC  1 (kslx,co)\j  a:xx((o)fi{xx!(ks2±,(o)  +  ayy(co)p{0y\ks21,co) 


&°Xk(s2x-S  1±),a») 

X - - - expti&lra-Tq)]. 

S(0)(0,o>) 


(3.11) 


The  assumption  that  the  correlation  coefficients 

u®Xp'  ,w)  ( i=x,y )  are  “fast”  functions  of  p'  across  the 

11  _ 

source  implies  that  their  Fourier  transforms  (f ' ,  co)  are 
“slow”  functions  of  f'.  We  can,  therefore,  make  the  ap¬ 
proximations 

pfJ(kslx,cj)  ~  p^(ks2x,to)^ffxJ(k(s2x  +  s1A)/2,w), 

(3.12a) 

gfyiksi  |  ,w)  «  gyy(ks2  ,,&))  =  g}yJ(k(s.2  I  +  Sj  ,  )/2,C>), 

(3.12b) 

and  it  follows  that  the  first  factor  in  Eq.  (3.11)  can  be  ap¬ 
proximated  by  unity.  We  then  obtain  for  the  complex  de¬ 
gree  of  coherence  of  the  far  field  the  formula 


Sm(M  s2x-s1±),a>) 

V "  (riSi,r2s2,ftj)  = - - - exp[i&(r2  -  /q)]. 


This  formula  expresses  another  reciprocity  relation 
pertaining  to  radiation  from  uniformly  polarized  planar 
QH  electromagnetic  sources,  namely,  the  spectral  degree 
of  coherence  ?7^(r1s1,r2s2,(u)  of  the  far  field  generated  by 
an  electromagnetic,  uniformly  polarized,  QH  source  is, 
apart  from  the  phase  factor  k(r2-rf),  proportional  to  the 
four-dimensional  spatial  Fourier  transform  of  the  spectral 
density  of  the  electric  field  in  the  source  plane .16 

We  can  summarize  this  part  of  the  analysis  by  saying 
that  we  have  shown  that  the  two  reciprocity  relations  as- 
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sociated  with  beams  produced  by  scalar  planar  QH 
sources  (Ref.  1,  Sec.  5.3.2)  have  analogs  for  beams  gener¬ 
ated  by  planar  uniformly  polarized  QH  electromagnetic 
sources. 


4.  POLARIZATION  PROPERTIES  OF  THE 
FAR  FIELD 

The  full  description  of  the  polarization  properties  of  a  ran¬ 
dom  electromagnetic  beam  at  some  point  r  requires  the 
specification  of  the  beam’s  degree  of  polarization  and  of 
the  state  of  polarization  of  its  polarized  portion. 

The  spectral  degree  of  polarization  of  a  random  electro¬ 
magnetic  beam  at  a  point  r=rs  in  the  far  zone  is  given  by 
the  expression 
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V4  Det  W("'Hrs,rs,  oj) 

1- - Z - (4-la) 

[Tr  W^"*(rs,rs,o>)]2 

or,  more  explicitly,  by  the  expression 
P^\rs,u>) 

yJ[W^\rs,rs,(o)  -  W^rs.rsM]2  +  4|W^)(rs,rs,<o)|2 

W£4(rs,rs,co)  +  W^\rs,rs,co) 

(4.1b) 

On  substituting  into  Eq.  (4.1b)  from  Eq.  (3.6)  for  the  ele¬ 
ments  of  the  matrix  W,  A)  of  the  far  field  generated  by  a 
uniformly  polarized  electromagnetic  QH  source  evaluated 
at  a  point  r1s1=r2s2  =  rs,  one  finds  that  the  spectral  de¬ 
gree  of  polarization  of  the  far  field  is  given  by  the  formula 


P(”Vs,&>)  = 


\j[axx(co)f^J(ks±,w)  -  ayy(w)fi'yy’(ks±,w)Y  +  4[«xy(«)^J(£s  , ,  w)] 
axx(to)f&(k8  x,o>)  +  ayy((o)/^}(k  S_L,  to) 


(4.2) 


Using  definition  (2.16),  formula  (4.2)  may  be  rewritten  in  the  form 


PM(rs,(o)  = 


V [v-xJ(ksL,(u)  -  S±M]2  +  4aM[/4°)(&S±M]2 

jti,{°J(ks±,co)  +  a(<o)p$(kB±,a>) 


(4.3) 


This  formula  expresses  the  spectral  degree  of  polarization 
of  the  far  field  generated  by  a  uniformly  polarized  planar 
QH  electromagnetic  source  in  terms  of  the  four¬ 
dimensional  spatial  Fourier  transforms  of  the  correlation 
coefficients  p^\p’  ,&>)  of  the  electric  field  in  the  source 
plane  and  the  factor  a(to),  defined  by  Eq.  (2.13a). 

The  normalized  Stokes  parameters  of  the  far  field 
s(°°'  (z  =  1,2,3)  are  related  to  the  elements  of  the  cross- 
spectral  density  matrix  W<00)  in  the  same  way  as  in  the 
source  plane.  For  uniformly  polarized  QH  sources,  appli¬ 
cation  of  Eqs.  (2.16)  and  (3.6)  leads  to  the  far-field  normal¬ 
ized  Stokes  parameters  as18 


fi^iks  |  ,(o)  -  (ks±,(o) 

jj.^(k s±,&))  +  a((o)jlfJ(ksL,(o)  ’ 


(4.4a) 


S2°Vs,<u)  : 


2\a(a)jRe[p^(ks  j  ,&>)] 

PtaO^rM  +  “MaC*  (&s±M  ’ 


(4.4b) 


2  V  a(w)lm\Jj${ksL,  &))] 

- 1 - .  (4.4c) 

Px °>(ks±,(o)  +  a(u>)p$(k  s±,&>) 

To  summarize,  formulas  (4.3)  and  (4.4)  express  the 
spectral  degree  of  polarization  and  the  normalized  spec¬ 
tral  Stokes  parameters  of  the  far  field  generated  by  a  pla¬ 
nar  uniformly  polarized  QH  source  in  terms  of  the  Fourier 
transforms  Zt<0),  u}°\  and  fi(0)  of  the  correlation  coeffi- 
cients  of  the  field  components  in  the  source  plane  and  in 
terms  of  the  ratio  a(oj)  of  the  spectral  densities  S[°\p,(o) 
and  S  ( p,(o )  of  the  components  of  the  electric  vector  in 
that  plane  [Eq.  (2.13a)]. 


SjTVs,  &>) 


5.  EXAMPLE:  THE  FAR  FIELD  GENERATED 
BY  A  MODEL  SOURCE 

We  will  illustrate  our  main  results  by  determining  the 
spectral  density,  the  spectral  degree  of  coherence,  the 
spectral  degree  of  polarization,  and  the  normalized  Stokes 
parameters  of  the  far  field  produced  by  sources  generated 
by  a  technique  described  in  Ref.  19  (see  also  Ref.  20).  Ac¬ 
cording  to  that  technique,  the  so-called  electromagnetic 
Gaussian  Schell-model  source  can  be  synthesized  from 
two  coherent,  linearly  polarized  plane  waves  transmitted 
through  two  mutually  correlated  phase-only  liquid-crystal 
spatial  light  modulators  (SLMs),  placed  in  the  arms  of  a 
Mach-Zehnder  interferometer. 

It  was  shown  in  Ref.  19  that  the  elements  of  the  cross- 
spectral  density  matrix  of  the  field  generated  by  such  a 
device  has  the  form  of  our  Eq.  (2.15),  with  the  spectral 
density  S(0'  and  the  correlation  coefficients  being 
Gaussian  functions,  viz., 


S(0\p,(o)  = 


AM 

—exp 


2cr2M 


0  =  x,y),  (5.1) 


P^M-PiM 


RyMexp 


IP2-P1I 

2<fM 


( i  =  X,y;j  =  x,y).  (5.2) 

Here  the  coefficients  Bjj(ai)  are,  in  general,  complex  num¬ 
bers.  Since  in  this  example  the  variances  <tx(co )  and  cry(co) 
are  the  same,  we  denoted  each  of  them  by  cr(co).  This  con¬ 
dition  ensures  that  the  source  is  uniformly  polarized.21 


2552  J.  Opt.  Soc.  Am.  A/Vol.  22,  No.  11 /November  2005 


Korotkova  et  al. 


There  are  various  constraints  on  the  parameters  enter¬ 
ing  expressions  (5.1)  and  (5.2).  First,  the  general  proper¬ 
ties  (e.g.,  nonnegative  definiteness)  of  the  cross-spectral 
density  matrix  of  the  field  in  the  source  plane  impose  cer¬ 
tain  conditions  on  the  parameters  of  the  source.  The  con¬ 
ditions  have  been  derived  in  Ref.  22  for  electromagnetic 
Gaussian  Schell-model  sources.  In  addition,  in  order  to 
ensure  that  the  field  produced  by  such  a  source  is  beam¬ 
like,  the  following  conditions  must  be  satisfied  (see  Ref. 
23): 


An2 

S(0)(f,  &))  =  - exp 

4-77 


(5.5) 


^2 

My0)(f',")  =  -Lb,j  exp 

LIT 


2 


(5.6) 


On  substituting  from  Eqs.  (5.5)  and  (5.6)  into  Eq.  (3.7) 
we  obtain  the  following  expression  for  the  spectral  density 
of  far  field: 


1  1  C2  2TT2 

4o2(a))  <5^.(cu)  2a ?  X2 


S(c°\rs,aj) 


Ak V 

2r2[l  +  a] 


l4x  exp(-  C^202/2) 


1  1  C2  2T72 

+  =  (5‘3) 

From  now  on  we  will  omit  the  explicit  dependence  of  all 
the  parameters  on  the  frequency.  In  order  that  the  elec¬ 
tromagnetic  Gaussian  Schell-model  source  generated  by 
means  of  the  system  of  the  SLMs  will  be  quasi- 
homogeneous, 

<r>  Sij  ( i=x,y ,  j  =  x,y).  (5.4) 

To  determine  the  behavior  of  the  far  field,  we  evaluate 
first  the  four-dimensional  Fourier  transform  of  the  spec¬ 
tral  density  (5.1)  and  of  the  correlation  coefficients  (5.2). 
One  readily  finds  that 


+  “4y  exP(-  (%yk2(P/2)\.  (5.7) 

As  before,  8,  the  angle  that  the  vector  r=rs  makes  with 
the  z  axis,  is  assumed  to  be  small,  so  that  cos  8~  1,  |s±|2 
=  sin2  8~  O2.24 

On  substituting  from  Eq.  (5.5)  into  formula  (3.13)  we 
obtain  for  the  spectral  degree  of  coherence  of  the  far  field 
the  expression 

97(”)(r1s1,r2s2,(o)  =  exp[-  a2^ |s2±  -  s1±|2/2]exp[t£(r2  -  rj)]. 

(5.8) 

From  Eqs.  (5.6)  and  (4.3)  we  readily  obtain  the  following 
expression  for  the  spectral  degree  of  polarization  of  the 
far  field: 


P(co)(rs,  to) 


exp(-  <?xxk2d2/2)  -  aS^y  exp(-  S^yk2  (P/2)]2  +  4:a^,y\Bxy\2  exp(-  S^yk2/)2) 
4*  exp(-  (%xk2(P/2)  +  ai?yy  exp(-  S^yk2(P/2) 


(5.9) 


Finally,  it  follows  from  Eqs.  (4.4)  and  (5.6)  that  the  nor¬ 
malized  spectral  Stokes  parameters  of  the  far  field  gener¬ 
ated  by  a  planar  electromagnetic  Gaussian  Schell-model 
source  are  given  by  the  expressions 

M  4c  exP(“  %xk2eP/2)  -  exp(-  tfyk2(P/2) 

1  ’  (%xexp(-  t%xk2(P/2)  +  ad^yexp(-  t%yk2(P/2)’ 

(5.10a) 


S  2°  Vs,  to) 


2  \  Q'4y|-®iylCOS  <Pxy  exp(-  (^.yhP/Pfe) 

<?xx  exp(-  J%xk2(P/2)  +  aS^y  exp(-  S^yk2(P/ 2)  ’ 

(5.10b) 


2  \  u 4y \Bxy | sin  (pxy  exp(-  4yk2^/2) 

S'  (rs  — - 

3  ’  4c  exP(_  4*V<92/2)  +  exP(_  <%yk2(P/2)  ’ 

(5.10c) 

where  (pxy  is  the  phase  of  the  coefficient  Bxy  [see  Eq.  (5.2)]. 

Figures  3-5  show  the  behavior  of  the  spectral  density, 
normalized  by  its  axial  value,  of  the  spectral  degree  of  po¬ 
larization  and  of  the  normalized  Stokes  parameters  of  the 
electric  field  in  the  far  zone,  generated  by  a  uniformly  po¬ 
larized  QH  electromagnetic  source,  specified  by  Eqs.  (5.1) 
and  (5.2). 


In  Figs.  3  and  4  the  off-diagonal  elements  of  the  cross- 
spectral  density  matrix  are  both  zero,  i.e.,  Bxy=Byx=  0.  We 
see  that  although  the  source  is  unpolarized,  the  far  field  is 
partially  polarized  and  the  degree  of  polarization  varies 
with  the  angle  0  and  also  with  the  difference  Sxx  -  Syy  be¬ 
tween  the  rms  widths  of  the  source  correlation  coefficients 
[see  Eq.  (5.2)].  One  can  show  from  Eq.  (5.9)  that  when 
Bxy=Byx= 0,  then  P(“)(rs,&))  =  |s<1^(rs,w)|,  while  S^Vs, cu) 
=  Sg“,(rs,&))  =  0.  Figure  3  pertains  to  the  case  for  which 
changes  of  the  degree  of  polarization  in  the  far  zone  occur 
only  near  the  beam  edge,  while  Fig.  4  illustrates  the  case 
for  which  changes  in  the  degree  of  polarization  occur 
within  the  beam  cross  section. 

Figure  5  illustrates  the  far-zone  properties  of  a  field 
generated  by  a  planar,  uniformly  polarized,  quasi- 
homogeneous  Gaussian  Schell-model  source,  which  is 
characterized  by  a  cross-spectral  density  matrix 
W^Vs  ,rs,  to)  whose  off-diagonal  elements  are  not  neces¬ 
sarily  zero.  Unlike  sources  with  “diagonal”  cross-spectral 
density  matrices  (Figs.  3  and  4),  the  field  generated  by 
this  source  cannot  be  treated  as  a  superposition  of  two  or¬ 
thogonally  polarized  scalar  fields.  It  is  evident  from  Eqs. 
(5.9)  and  (5.10)  that  the  state  of  polarization  of  the  far 
field  generated  by  a  QH  source  represented  by  a  cross- 
spectral  density  matrix  with  nonvanishing  offdiagonal  el¬ 
ements  has  nontrivial  dependence  on  angle  8. 
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The  spectral  density  and  the  polarization  properties  of 
the  far  field  generated  by  a  uniformly  polarized  QH 
source,  characterized  by  Eq.  (5.1M5.4),  possess  rotational 
symmetry.  Figure  6  illustrates  this  property  for  the  case 
represented  in  Fig.  5.  The  black-to-white  (gray-scale) 
background  in  this  diagram  corresponds  to  the  values  of 
the  spectral  density  of  the  far  field  (normalized  by  its 
axial  value)  in  the  range  from  0  to  1.  The  black  ellipses 
superimposed  on  this  background  represent  the  spectral 
polarization  ellipses  of  the  fully  polarized  portion  of  the 
far  field.  The  parameters  of  each  of  the  spectral  polariza¬ 
tion  ellipses  can  be  readily  calculated  from  the  Stokes  pa¬ 
rameters  [see  Eq.  (5.10)]  with  the  help  of  well-known  re¬ 
lations  (see  Ref.  13,  Sec.  10.9.3,  see  also  Ref.  25). 


APPENDIX  A:  NECESSARY  AND 
SUFFICIENT  CONDITIONS  FOR  A 
UNIFORMLY  POLARIZED  PLANAR  QH 
SOURCE 

We  will  show  that  a  QH  source  is  uniformly  polarized  at 
frequency  a>  if  and  only  if  the  spectral  densities  S{x\p,co) 
and  Sy  (p, oj)  are  proportional  to  each  other,  i.e.,  if 

S(y0\p,co)  =  a(co)S{x°\p,co)  (Ala) 

and  if,  in  addition,  the  correlation  coefficient  p°y(p,p,io) 
is  independent  of  position,  i.e.,  if 


Fig.  3.  The  spectral  density  S(rs,<i>)  and  the  spectral  degree  of  polarization  of  the  far  field  generated  by  transmission  of  a  linearly 
polarized  beam  through  a  system  of  spatial  light  modulators  that  produces  a  uniformly  polarized  QH  field  with  parameters  a=  1,  <x 
=  1  mm,  Bxy= 0,  Sxx  =  0.1  mm,  8yy=0.5  mm,  X  =  0.6328  ^tm.  The  spectral  density  is  normalized  by  its  on-axis  value. 


Fig.  4.  The  curves  represent  the  same  quantities  as  in  Fig.  3,  but  the  value  of  parameter  a  is  0.05. 
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Fig.  5.  (a)  The  degree  of  ploarization  P(rs ,  u>)  and  the  spectral  density  S(rs ,  a>)  of  the  far  field  generated  by  the  source  with  parameters 
a=0.44,  <r=  1  mm,  Bxy=0.525,  i5a.a.=0.15  mm,  Syy=  0.25  mm,  <5^=0. 3  mm,  \  =  0.6328  /im,  <pxy=\/6.  (b)  Normalized  Stokes  parameters  of  the 
far  field  generated  by  the  same  source  as  in  (a). 


p*?(p>p>")  =  /4y  (")• 


= 


(Alb) 


First,  if  conditions  (Ala)  and  (Alb)  hold,  then  one  finds, 
with  the  help  of  the  Eqs.  (2.10),  that  the  normalized 
Stokes  parameters,  given  by  the  expressions 


sfM  = 


1  -  a(w) 
1  +  a(co) 


(A2a) 


: 


S^M  = 


2VnMIm[/4yM] 

1  +  a(co) 


(A2c) 


are  independent  of  position.  Hence,  such  a  source  is  uni¬ 
formly  polarized. 

Conversely,  suppose  that  the  state  of  polarization  is 
uniform  across  the  source.  The  normalized  Stokes  param¬ 
eters,  given  by  the  formulas  (2.12),  are  then  independent 
of  p,  say, 


1  +  a((o ) 


(A2b) 


((  =  1,2,3) 


(A3) 


Korotkova  et  al. 


Vol.  22,  No.  ll/November  2005/J.  Opt.  Soc.  Am.  A  2555 


Equations  Eqs.  (2.12)  and  (A3)  then  imply  that 


S?\p,co)  l  +  sfM 
Sf(p,a>ri~S<0)M  = 


(A4) 


APPENDIX  B:  DERIVATION  OF  EQ.  (3.3) 

In  this  appendix  we  will  evaluate  the  four-dimensional 
Fourier  transform  of  the  cross-spectral  density  matrix  in 
the  source  plane,  viz., 


and  that 


Re[^(p,P,")] 


s!j0,M 

\/l-[s  i0,(o>)]2’ 


(A5a) 


Wp (pi ,  p2 ,  <u) exp[ -  t(fj  •  Pl 


+  f2  •  p2)]d2Pid2p2  (i  =  x,y,  j  =  x,y), 

(Bl) 


Mpx^p.p,")] 


sj,0)M 

vl-[sfV)]2' 


(A5b) 


From  Eq.  (A4)  it  follows  that  the  spectra  S'"  and  S'"  are 
proportional  to  each  other,  and  from  Eqs.  (A5)  it  follows 
that  the  correlation  coefficient  p}"v)(p1,p2,io),  evaluated  for 
Pi=p2  =  p,  is  independent  of  p.  Formulas  (A4)  and  (A5) 
also  show  how  the  ratio  S(p{p,to)IS'"\p,oj)  of  the  spectra 
and  the  correlation  coefficient  p  (p ,  p ,  cn)  can  be  deter¬ 
mined  from  the  knowledge  of  the  normalized  Stokes  pa¬ 
rameters. 


where 


W$\p1,p2,<o)  = 


Uy(«j)S(0i((pi  +  p2)/ 2,  (o)p.p(p2  -  p1;  io)  . 


(B2) 


Introducing  the  variables 


Ap  =  p2-pi,  p  =  (pi  +  p2)/2,  (B3) 


we  can  express  formula  (Bl)  in  the  form 


Fig.  6.  Illustration  of  rotational  symmetry  of  the  spectral  density  and  of  polarization  of  the  far  field  generated  by  the  same  source  as  in 
Fig.  5(a),  plotted  as  functions  of  polar  angle  6  in  the  interval  -3  X  10-3  rads  (?=£  3  X  10-3  rad.  The  black-to-white  (gray-scale)  background 
corresponds  to  the  normalized  values  between  0  and  1  of  the  spectral  density  of  the  far  field,  normalized  by  its  axial  value.  The  ellipses 
superimposed  on  this  background  are  the  spectral  polarization  ellipses  of  the  fully  polarized  portion  of  the  far  field.  Their  centers  are 
chosen  to  be  located  along  the  two  mutually  orthogonal  directions  at  points  where  the  spectral  density  normalized  by  its  axial  value  takes 
on  the  values  1,  0.85,  0.6,  0.4,  0.25,  and  0.1. 
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“y(") 


")/4;0)(Ap,u>)exp[-  i(fi  •  (p  -  Ap/2)  +  f2  •  (p  +  Ap/2))]d2pd2(Ap) 


Wff)(f1,f2>«)  =  ^  j  j  S<°>(p, 

■  j  j  S(0>(p, o»)49)(A p, 0>)exp[-  i(p  ■  (fj  +  fa)  +  Ap  •  (f2  -  f1)/2)]d2pd2(Ap) 


(2tt) 


2  J  S(0Hp,v)exp[-  ip  (f1  +  f2)]d2p-^—^  J  p[;0)(Ap,(o)exp[-  iAp  ■  (f2  -  fj)/2]d2(Ap) 


:  +  f2, o,)4?>((f2  -  fJ/2, oi) . 


(B4) 


This  formula  is  analogous  to  the  corresponding  expres¬ 
sion  that  pertains  to  the  scalar  case  [see  Eq.  (5.3-19)  in 
Ref.  1], 
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